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The ACDM cosmological model is a well defined, simple and predictive model which is consistent 
with the majority of current cosmological observations. Despite of these successes there are specific 
cosmological observations which differ from the predictions of ACDM at a level of 2a or higher. These 
observations include the following: 1. Large Scale Velocity Flows (ACDM predicts significantly 
smaller amplitude and scale of flows than what observations indicate), 2. Brightness of Type la 
Supernovae (Snia) at High Redshift z (ACDM predicts fainter Snia at High z), 3. Emptiness of 
Voids (ACDM predicts more dwarf or irregular galaxies in voids than observed), 4. Proflles of Cluster 
Haloes (ACDM predicts shallow low concentration and density profiles in contrast to observations 
which indicate denser high concentration cluster haloes) 5. Profiles of Galaxy Haloes (ACDM 
predicts halo mass proflles with cuspy cores and low outer density while lensing and dynamical 
observations indicate a central core of constant density and a flattish high dark mass density outer 
proflle), 6. Sizable Population of Disk Galaxies (ACDM predicts a smaller fraction of disk galaxies 
due to recent mergers expected to disrupt cold rotationally supported disks). Even though the origin 
of some of the above challenges may be astrophysical or related to dark matter properties, it should 
be stressed that even on galactic and cluster scales, the effects of dark energy on the equilibrium and 
stability of astrophysical systems are not negligible and they may play a key role in the resolution 
of the above puzzles. Here, I briefly review these six challenges of ACDM and discuss the possible 
dark energy properties required for their resolution. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



1. INTRODUCTION 



Accumulating diverse observational evidence have indicated that the universe has entered a phase of accelerating 
expansion. Such observations include direct geometrical probes (standard candles like Snia P, H, H, 01 , gamma ray 
bursts and standard rulers like the CMB sound horizon [^Q) and dynamical probes (growth rate of cosmological 
perturbations JSj probed by the redshift distortion factor or by weak lensing ^]). 

All these observational probes are converging towards confirming the accelerating expansion of the universe assuming 
the homogeneity of the universe. They have ruled out at several a a flat matter dominated universe and they have 
produced excellent fits for the simplest cosmological model predicting accelerating cosmic expansion. This model is 
based on the assumptions of flatness, validity of general relativity, the presence of the cosmological constant A and 
Cold Dark Matter (ACDM 

From the theoretical viewpoint the main weak points of ACDM include p^ : 



The Fine Tuning Problem: What is the physical mechanism that sets the value of A to its observed value which 
is 120 orders of magnitude smaller than the physically anticipated value? 

The Coincidence Problem: Why is the energy density corresponding to the cosmological constant just starting 
to dominate the universe at the present cosmological time? 



Despite of efforts to increase the complexity of ACDM (using eg quintessence [ll| or modified gravity[l3|) in order to 
address the above weak points there has been no successful alternative that addresses the above problems without 
replacing them with other similar ones involving fined tuned parameters. Since the theoretical weaknesses of the 
model have lead to no successful alternative it may be useful to identify the observational weak points of ACDM and 
use these as a guide to building alternative models. 

In view of the fact that ACDM is a simple, well defined and predictive model, it is important and straightforward to 
test its validity using a wide range of observational probes. If some of these observational probes indicate inconsistency 
of ACDM with observations then it is interesting to consider the modifications of the model required to establish 
consistency with observations. 

Most approaches in testing the consistency of ACDM with observations have focused on comparing ACDM with 
alternative models or parameterizations on the basis of a bayesian analysis using the geometrical and dynamical 
probes mentioned above [EH, B 3, Hi B, B ■ Due to its simphcity and acceptable quality of fit, ACDM usually 
comes out as a winner in such a comparison \13^ . 

Despite of the simplicity and apparent consistency of ACDM with most cosmological observations there are specific 
observational challenges for the model which have developed and persisted during the past few years. Some of 
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these challenges involve galactic scale phenomena and it has been common wisdom that they will be resolved once 
astrophysical effects on these scales are better understood. Other challenges however, involve phenomena on scales 
larger than ~ 10h~^Mpc and these may require more drastic modifications of the model in order to be resolved. 
Such large scale challenges of ACDM include the observed high amplitude of large scale velocity flows on scales 

- woh-^Mpc [11, lillii, the unexpected brightness of high redshift Type la supernovae (Snla)[l3|, the halos of 
massive clusters of galaxies which are more concentrated and denser than predicted by ACDM [l^ and the emptiness 
of voids which is unexpected in the context of ACDM On smaller (galactic) scales ACDM is challenged 



by observations of constant density galactic halo cores instead of the ACDM predicted cuspy central cores [24I , the 
higher than expected density of outer galactic haloes 23] and the sizable population of cold rotationally supported 
disk galaxies [2J|. 

Since the above effects are statistically significant at 2a level or more it is unlikely that they are all statistical 
fluctuations. In fact, it is possible that the resolution of the above puzzles will require more than a better understanding 
of astrophysical effects present on galactic scales. It may require a significant modification of the cosmological scale 
properties of the standard ACDM model such as the properties of gravity, dark energy or dark matter. 

The goal of the present paper is to review the above phenomena challenging the foundations of the standard ACDM 
cosmological model. I will also discuss possible features of the model that may require modification in order to improve 
consistency with the above observations. 

It should be stressed that this is not a complete list of cosmological puzzles related to the standard ACDM cos- 
mological model. There are other challenges related to the statistical isotropy of the CMB and the Axis of Evil [25j 
(anomalous alignment of CMB multipoles in the direction I ~ —100°, b — 60°) which may be less related to the 
properties of dark matter or dark energy. Such challenges are not discussed in the present brief review even though 
they may be related to the high amplitude and coherence bulk flows discussed in the next section. 



2. CHALLENGING ACDM 
A. Large Scale Velocity Flows 

The bulk flow corresponding to the CMB dipole is closely related to the amplitude and growth rate of fluctuations 
on large scales, and can be used to test cosmological models [1^ . A number of large scale velocity surveys have been 
undertaken [27| in the past two decades and a significant amount of peculiar velocity data on a wide range of scales is 
currently available. The issue of comparing such sparse surveys with expectations from cosmological models has also 
been investigated by several studies |28| . 

A combined sample of peculiar velocity data has been recently used [3, [l^ to investigate the amplitude and 
coherence scale of the dipole bulk flow. It was found that the dipole moment (bulk flow) of the combined sample 
extends on scales up to lOOh^^AIpc (z < 0.03) and perhaps up to QOOh^^Mpc {z < 0.2 [15|]) with amphtude 
larger than 400/cm/sec (perhaps up to lOOOkm/sec [lal)- The direction of the flow has been found consistently 
to be approximately in the direction I ~ 285°, b ~ 10° in rough agreement with the CMB dipole {I ~ 276°, b ~ 30°). 
Similar results implying a large bulk low ~ 500km/ sec on scales up to lOOh^^Mpc were recently also obtained in 
Ref. dll. The expected rms bulk flow in the context of ACDM normafized with WMAP5 (f^om, ^s) = (0.258, 0.796) 
on scales larger than 50h~^Mpc is approximately 110km/ sec while the probability that a flow magnitude larger than 
AOOkm / sec is realized in the context of the above ACDM normalization on scales larger than 5Qh~^AIpc is less than 
1%. 

This is also demonstrated in Fig. 1 (from Ref. [IJ]) where the (fiom, cs) confidence contours obtained from the 
observed velocity flows (dashed lines) are superposed with the corresponding contours obtained from WMAP5 data 
(blue solid fines) and from WMAP5+Baryon Acoustic OsciUations+SnIa (WMAP5+BA0+SN: red dashed fine). The 
probability of consistency of bulk flow data with ACDM would be even lower if the data of Ref. [lB| were considered 
where a flow of more than 600km/ sec was observed on scales of 600h~^Mpc. 

A potential resolution of the above described conflict between the high z WMAP5 normalization of ACDM and the 
low z normalization implied by the observed bulk flows could involve the existence of superhorizon sized non-Gaussian 
and non- inflationary inhomogeneities [1^ , a large void at distances of order gigaparsecs ^3] , or a redshift dependent 
(T8 which changes by a factor of 2 between high z and low z due to an unknown physical reason. Other possibilities 
include a very large statistical fluctuation, a redshift dependence of Newton's constant or a redshift dependence of the 
dark energy equation of state parameter w — w{z) leading to amplified gravity and dark energy clustering at early 
times (u>(z) > — 1 at 2; > 0.2). 
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FIG. 1: The (nom,o"8) confidence contours obtained from tfie observed velocity flows (dashed lines) are superposed 
with the corresponding contours obtained from WMAP5 data (blue solid lines) and from WMAP5+BAO+SN (red dashed line) 
(fromRef. Q). 
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FIG. 2: The UnionOSQ distance moduli data superposed with the best fit ACDM model (Qom ~ 0.29) dashed line and with 
the best fit {wo,wi) = ( — 1.4,2) model (f2om = 0.30) continous line. Notice that at high redshifts z the distance moduli tend 
to be below the ACDM best fit while the trend is milder in the PDL crossing best fit model (from Ref. [l7l |. 



B. Bright High z Snia 

As discussed in the introduction, geometrical tests of ACDM usually involve a bayesian comparison of ACDM with 
other dark energy parametrizations. This approach has not revealed so far any statistically significant weak points of 
the model with respect to the geometrical and dynamical probes considered. 

Apart from the bayesian analysis approach, the ACDM model can be tested by comparing the real SnIa data with 
Monte Carlo simulations consisting of fictitious cosmological data that would have been obtained in the context of a 
ACDM cosmology. This comparison can be made on the basis of various statistics which attempt to pick up features 
of the data that can be reproduced with difficulty by a ACDM cosmology [l7'|. The existence of such features is hinted 
by the form of the likelihood contours in various parameter planes containing parameter values corresponding to flat 
ACDM. For example, most SnIa datasets producing likelihood contours in the Q\ — parameter plane have the 
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1(7 contour barely intersect the line of flatness J^a + = 1 at the lower left side of the contour [3, Q]- Similarly, 
likelihood contours based on either Snia standard candles or standard rulers (CMB sound horizon or Baryon Acoustic 
Oscillations (BAO)) and constraining the parametrization [sH 



w{z) — Wo + wi 



1 + z 



(2.1) 



systematically have the point corresponding to ACDM (wo ) — (—1,0) at the lower right edge of the la contour 
while the best fit involves wo<— 1, '*i'i>o[l,|3, Is^ IbsI . This feature has persisted consistently over the last 
decade and over different accelerating expansion probes [3^ (SnIa standard candles and CMB-B AO standard rulers) . 
Even though the statistical significance of these features when viewed individually is relatively low, their persistent 
appearance makes it likely that there are systematic differences between the cosmological data and ACDM predictions. 



Tested Mo(lel:w(]=- 1 .4, v 




N,^ 105 from Gold06 data (z,.=0.47) 




FIG. 3: a: A histogram of the probabihty distribution of Nmc obtained using Monte Carlo ACDM data (^lom = 0.34) in the 
context of the Gold06@| dataset. The thick green dashed line corresponds to the crossing redshift Zc of the real GoldOG data, 
b: Similar histogram for the PDL crossing model (wojWi) = ( — 1.4,2) (best fit ilom = 0.34) instead of ACDM . Notice that the 
crossing redshift Zc corresponding to the real Gold06 data is a much more probable event in the context of this cosmological 
model (from Ref.T^). 

One such difference in the context of SnIa data has been recently pointed out by Kowalsky et. al. [3| where it was 
stated that there is 'an unexpected brightness of SnIa data at z > 1'. This feature is even directly visible by observing 
the SnIa distance moduli superposed with the best fit ACDM model (dashed line in Fig. 2) where most high z moduli 
are below the best fit ACDM curve (obviously the reverse happens at low redshifts to achieve a good fit) . Notice that 
this bias is smaller in the context of a parametrization that crosses the PDL w — —1 (continuous line in Fig. 2). ^ 

This anomalous behavior of the data with respect to the ACDM best fit may be attributed to the systematic 
brightness trend of high redshift SnIa with respect to the best fit ACDM model. It is likely that this bias of the SnIa 
data with respect to ACDM best fit is also responsible for the systematic mild preference (at la) of the SnIa data for 
a w(z) crossing the w — —1 line. 

In order to study quantitatively the likelihood of the existence of the above described bias in the context of a 
ACDM cosmology, we may use a statistic flT'l (the Binned Normalized Differences (BND)) specially designed to pick 
up systematic brightness trends of the SnIa data with respect to a best fit cosmological model at high redshift. The 
BND statistic is based on binning the normalized differences between the SnIa distance moduli and the corresponding 
best fit values in the context of a specific cosmological model (eg ACDM). These differences are normalized by 
the standard errors of the observed distance moduli. We then focus on the highest redshift bin and extend its size 
towards lower redshifts until the Binned Normalized Difference (BND) changes sign (crosses 0) at a redshift Zc (bin 
size Nc). The bin size Nc of this crossing (the statistical variable) is then compared with the corresponding crossing 
bin size Nmc for Monte Carlo data realizations based on the best fit model. It may be shown [TtI that the crossing 
bin size Nc obtained from the UnionOS and GoldOG data with respect to the best fit ACDM model is anomalously 
large compared to N^c of the corresponding Monte Carlo datasets obtained from the best fit ACDM in each case. 
In particular, only 2.2% of the Monte Carlo ACDM datasets are consistent with the Gold06 value of Nc (see Fig. 



^ In the PDL crossing model wc fix wo, wi and vary Qom only, in order to mimic the ACDM number of parameters. 
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3a) while the corresponding probabihty for the UnionOS value of Nc is 5.3%. Thus, according to this statistic, the 
probability that the high redshift brightness bias of the UnionOS and GoldOG datasets is realized in the context of 
a {wo,wi) = (—1,0) model (ACDM cosmology) is less than 6%. The corresponding realization probability in the 
context of a {wq, wi) — (—1.4, 2) model is more than 30% for both the UnionOS and the GoldOG (see Fig. 3b) datasets 
indicating a much better consistency for this model with respect to the BND statistic. 

This result reveals a potential challenge for ACDM cosmology and provides the motivation for obtaining additional 
Snia data at high redshifts z > 1 which may confirm or disprove the anomalous high z Snia brightness which is mainly 
responsible for the low probability of the high z Snia data in the context of KCDM . 

Clearly, the unexpected high z brightness of Snia can be interpreted either as a trend towards more deceleration at 
high z than expected in the context of ACDM or as a statistical fluctuation or finally as a systematic effect perhaps 
due to a mild Snia evolution at high z. However, in view of the fact that a similar mild trend for more deceleration 
than expected at high z is also observed in the context of standard rulers [1, [H, [s^ , the latter two interpretations 
are less likely than the first. 



C. The Emptiness of Voids 

Cosmological simulatio ns p erformed in the context of ACDM predict [13] that many small dark matter haloes 
should reside in voids [Tol l20l [2^1]. This is consistent with observations on large scales involving giant voids defined 
by IO^^Mq haloes [s^. Smaller voids however {r^ lOMpc) look very empty. Dwarf galaxies do not show a tendency 
to fill these voids even though ACDM predicts that many dwarf dark matter haloes should be in the voids. In fact 
as discussed in [2l| , if efhciency of conversion of ACDM halos to galaxies observable in optical or HI emission were 
independent of environment then we would expect that about ten galaxies with — IS < Afs < — 10 (extreme dwarfs) 
are in the Local Void while none is observed [36j. 

This dwarf galaxy overabundance problem of ACDM (the 'void phenomeno n ' 13711 ) is also connected to the predicted 
number of dwarf satellites in the Local Group (the 'missing satellite problem'jS^): the theory predicts a factor often 
more haloes as compared with the observed number of dwarf galaxies. For example the ACDM model predicts that 
thousands of dwarf dark matter haloes should exist in the Local Group [H, [ll,!!^!, while only ^ 50 are observed. 
Recent discoveries of very low luminosity dwarfs [JTl and careful analysis of incompleteness effects in SDSS [U [Z^l 
bring the theory and observations a bit closer, but the mismatch seems is still present. 

Potential resolutions of the above tension between ACDM theory and observations involve incompleteness of obser- 
vational sample, failure of many dwarf haloes to form stars[43] in the context of a mass dependent bias model (Halo 
Occupation Distribution [44]) or peculiar properties of dark matter and/or dark energy which accelerate growth of 
perturbations and allow gravity to clean up voids at early times. Even though the simplest resolution of the 'void 
phenomenon' could involve the use of a simple bias model where galaxy formation is driven predominantly by the 
mass of the host dark matter halo [i^ , the effectiveness of this approach at the level of dwarf galaxy luminosities is 
under debate [Tol. [2H and it is possible that an environment dependent bias may be required. 



D. Galaxy Halo Profiles 

The ACDM theory predicts that dark matter halos have a specific density distribution that follows the well-known 
Navarro, Frenk, White (NFW) [H, ^ profile: 

Pnfw{R) = j-^i — w-.^", p, ^2 (2-2) 

where and ps are the characteristic radius and density of the distribution. A useful parameter characterizing the 
profile is the concentration parameter c defined as is c = r^ir/r^ where ryir is the virial radius of the system, and 
Ps are related to each other (e.g. 113]), so eq. (12. 2p is rather a one-parameter family of profiles. 

A quite remarkable number of observations show that NFW profiles, displaying an inner cusp, are inconsistent with 
data[4a]. In fact, the latter indicate profiles with a different characteristic, a central density core, i.e. a region where 
the dark matter density remains approximately constant. 

In addition to the above well-known evidence for which in the inner regions of galaxies (i? < 2rd where is the 
stellar disk radius) the dark matter haloes show a flatfish density profile, with amplitudes up to one order of magnitude 
lower than the ACDM predictions, at outer radii {R > ird) the measured dark matter halo densities are found higher 
than the corresponding ACDM ones. The dark matter halo density, known to have a core in the internal regions, does 
not seem to converge to the NFW profile at 4 — 6rd [1^ . This implies an issue for ACDM that should be investigated 
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in the future, when, due to improved observational techniques, the kinematic information will be extended to the 
lOOfcpc scale. 

A possible resolution of the puzzle of higher than expected dark matter halo density in the galactic haloes is that 
massive halos themselves were assembled at high redshift [20*1 . If this is the case, modifying the properties of dark 
energy could play a role in shifting the epoch of galaxy formation towards earlier times. Alternatively, modified 
gravity theories or clustering of dark energy may also be considered as a potential resolution of this puzzle. 

E. Cluster Halo Profiles 

In the ACDM context, detailed N-body simulations have established a clear prediction that CDM-dominated 
cluster halos should have relatively shallow, low-concentration mass profiles, where the logarithmic gradient flat- 
tens continuously toward the center with a central slope tending towards r~^, interior to a characteristic radius, 

r, ~ 100 - 2ookpc -h-^ [H, EE H mi, [iljp . 

Multiply-lensed images of various clusters [18|| have been used to derive the inner mass profile [l^ljWith the outer 
profile determined from weak lensing [HH]. Together, the full profile has the predicted NFW form [4^, but with a 
surprisingly high concentration c = and high density when compared to the shallow profiles of the standard 
ACDM model [s^ [sH] • This result is verified by using not only the lensing based mass profile but also the X-ray and 
dynamical structure in model independent analyses (57j . 

A potential resolution of the above discrepancy between observed cluster profiles and ACDM predictions is that the 
central region of clusters collapsed, as in the case of galaxies, earlier than expected ie at z > I, significantly earlier 
than in the standard ACDM , for which clusters form at z < 0.5. 

The presence of massive clusters at high redshift ~ 2), and the old ages of their member galaxies [H, [59[, may 
also imply clusters collapsed at relatively early times 601, for which accelerated growth factors have been proposed, 
adopting a generalized equation of state for dark energy (61j . Such an equation of state would allow for a non- negligible 
dark energy density at early times. Thus, as in the case of galaxy formation, the properties of dark matter and/or 
dark energy could also play a significant role in the resolution of this puzzle. 

F. Overpopulation of Disk Galaxies 

Roughly 70% of Milky- Way size dark matter halos are believed to host late-type, disk dominated galaxies [63 |. 
Conventional wisdom dictates that disk galaxies result from fairly quiescent formation histories, and this has raised 
concerns about disk formation within the hierarchical ACDM cosmology [6^ l63l. Recent evidence for the existence 
of a sizeable population of cold, rotationally supported disk galaxies at z ~ 1.6 [65] is particularly striking, given that 
the fraction of galaxies with recent mergers is expected to be significantly higher at that time ,66.] . High-resolution, 
dissipationless N-body simulations [63] studying the response of stellar Milky- Way type disks to such common mergers 
show that thin disks do not survive the bombardment. The remnant galaxies are roughly three times as thick and 
twice as kinematically hot as the observed thin disk of the Milky Way. However, despite of such indications a real 
evaluation of the severity of the problem is limited by both theoretical and observational concerns. 

The role of dark energy in the resolution of this and other astrophysical scale puzzles should not be underestimated. 
For example, it has been demonstrated that the effects of dark energy on the equilibrium and stability of astrophysical 
structures is not negligible, and can be of relevance to describe features of astrophysical systems such as globular 
clusters, galaxy clusters or even galaxies [gi, H, H [nl • It has recently been demonstrated that the dark energy fluid 
changes certain aspects of astrophysical hydrostatic equilibrium. For example, the instability of previously viable 
astrophysical systems when dark energy is included has been demonstrated as due to the repulsive non local dark 
energy force acting on the matter distribution [72]. With the proper evolution of the dark energy equation of state, 
this repulsive force may also lead to a modification of the profile of the virialized structures thus addressing some of 
the above discussed puzzles on galactic and cluster scales. 

3. DISCUSSION - CONCLUSION 

I have reviewed six of the potential observational challenges for the ACDM cosmological model (as normalized by 
WMAP5) pointing out that there are such challenges on both large and small cosmological scales. The observations 
conflicting the WMAP5 normalized ACDM model at a level of 2a or larger include the following: 

• Large Scale Velocity Flows (ACDM predicts significantly smaller amplitude and scale of flows than what obser- 
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vations indicate), 

• Brightness of Type la Supernovae (Snia) at High Redshift z (ACDM predicts fainter Snia at High z), 

• Emptiness of Voids (ACDM predicts more dwarf or irregular galaxies in voids than observed) , 

• Profiles of Cluster Haloes (ACDM predicts shallow low concentration and density profiles in contrast to obser- 
vations which indicate denser high concentration cluster haloes) 

• Profiles of Galaxy Haloes (ACDM predicts halo mass profiles with cuspy cores and low outer density while 
lensing and dynamical observations indicate a central core of constant density and a flatfish high dark mass 
density outer profile), 

• Sizable Population of Disk Galaxies (ACDM predicts a smaller fraction of disk galaxies due to recent mergers 
expected to disrupt cold rotationally supported disks). 

Even though some of the puzzles discussed here may be resolved by more complete observations or astrophysical 
effects, the possible requirement of more fundamental modifications of the ACDM model remains valid. 

It is interesting to attempt to identify universal features which connect these puzzles and could therefore provide 
a guide for their simultaneous resolution. The large scale coherent velocity flows along with the high density dark 
matter haloes for both galaxies and clusters seem to hint towards a more effective mechanism for structure formation 
at early times (z > 1) than implied by ACDM . This improved effectiveness could possibly be provided by a mild 
evolution of Newton's constant G (higher C? at z > 0.5) or by an evolution of the dark energy equation of state w such 

that w{z) > — 1 at z ~ 0.5 |61||. Both of these effects are expected to amplify structure formation at early times and it 
would be interesting to analyze quantitatively the predictions implied by the evolution of G or w with respect to the 
velocity flow and high dark matter density puzzles. The Bright High z Snia puzzle would also benefit significantly by 
a mild evolution of w or G which would imply stronger deceleration at z > 1 than implied by ACDM . 

The improved efficiency of gravity at early times could also help emptying the voids from dark matter haloes and 
their corresponding galaxies thus making theoretical predictions more consistent with observations. On the other 
hand, the increased gravitational acceleration would also produce higher peculiar velocities that could lead to more 
mass inside the voids. Therefore, the predicted emptiness of voids in models with an evolving G ox w requires a 
detailed study. 

In conclusion, the six puzzles for ACDM discussed in the present study provide a fertile ground for the development 
of both new theoretical model predictions on the corresponding observables and new observational data that would 
either establish or disprove these challenges for ACDM . 
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